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Abstract 

We present an algorithm for the construction of the branching functions in the vacuum 
sector for affine Lie algebras based on the string hypothesis solution to a system of Bethe 
equations for generalized RSOS models. We also mention how the ground state structure 
and features of the excitation spectra like the Brillouin zone schemes of these models (and 
those in the same universality classes) can be extracted from combinatoric arguments and 
encoded in Lie algebraic terms. 
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The exact solution of an exactly solvable model in 1+1 or 2+0 dimensions usually 
implies the parametrization of the dynamical quantities of interest (like the energy and 
momentum eigenvalues, and often, even the eigenvectors) in terms of the solutions of a 
coupled set of transcendental equations, called Bethe equations 0. In the critical case, 
they are usually written in terms of trigonometric functions, and for the generalized RSOS 
models, they are of the form 



sinh(^(Af + i(su; p \a a ))) ] N = ^ ^ Q sinh(^(A^ - Aj b) + i(a a \a b ))) 



smh(^ L (X ( - a) - i(su p \a a ))) 



N^NslC- 1 ], (2) 



where iV is the size of the lattice a a are the simple roots, uj p are the fundamental weights, 
(•|-) is the canonical bilinear form on the dual to the Cartan subalgebra, C is the Cartan 
matrix, L = I + g, for integer /, g is the dual Coxeter number and s characterizes the type 
of fusion. It should be mentioned that these have been derived in [[|only for the JKMO 
[jTjmodels that generalize the ABF ||models to the case (and their fusion hierarchies), 
and have only been postulated for the other cases. (See ||for a discussion.) Here is 
some phase factor that influences what the allowed solutions are. Since all the quantities of 
physical interest are parametrized by solutions to these equations (which are very difficult 
to solve), the standard practice is to assume a particular form for the roots called the 
"string hypothesis" [jljand hope that the equations that are set up in the thermodynamic 
limit to compute the spectrum or the free energy are consistent with this ansatz. In fact, 
one does more - one makes the tacit assumption that important physical quantities are 
"diagonalized in the string basis," and therefore classifying the spectrum in terms of the 
string form of the roots of the Bethe equations becomes significant. In what follows, the 
string hypothesis (as in @0 and 0) can be written as: 

x\ a) = A + it-\ 3 + 1 - 2k), 1 < h < J (3) 

with t~ x = (a a \a a )/2. X is real and denotes the center of the string and j denotes 
its length, which is further assumed to satisfy 1 < j ' < t a l. We shall therefore impose 
Y^j^i jMj = N a , where M- denotes the number of strings of colour (a), length j and 
iV a is given by (0). These assumptions were made in the bulk in 0(and subsequently in 
0and 0), but here we shall impose these on finite-size lattices. (Of course, the same 
thermodynamic limit can be recovered.) Despite the well-known examples of non-string 
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like behaviour (see ||, [lC]for a list of references) the results presented here do seem rather 
striking, and it would be interesting to understand to what extent these results are related 
to the actual models one would like to solve. 

As always, we multiply out the Bethe equations for the components of each string, 
and end up with equations for the real parts of the roots. We then take the logarithm of 
(Pso that the integer branches are distinct. For the simply-laced algebras, the formulas 
look simpler, (with all t a = 1) and for this case, we have 

l M (h) 

Nt^ s (X^) = 2*il$ + E ®% b) (K (a) - X t (b) ), (4) 

b=l k=l v=\ 

where Xjj a ^ labels the center of the (i-th string of length j and color (a) which is a root of 
eq. ([!]). The functions and are defined below. 

min(j,s) 

t^ s (X) = 5 ap /(A;|i-a|+2fc-l), (5) 

fe=i 

tmin(j,k) — l 
/(A; \j-k\)+2 <T /(A; \j - k\ + 2i) + f(X;j + k) 

tn(j,k) 



I ab [ £ /(A; \j - k\ + 2i - 1) 



i=i 



1 /sinh \ -r{in — A) 



/(A;n)= —In , - U , ■ ' , (7) 



2%i Vsinh i^(m + A) . 

for integer values of n/L, and is otherwise. I a b is the incidence matrix of the respective 
Dynkin diagrams. (For the non-simply laced cases, the branches are chosen following the 
same principle ||, but they look a little more complicated, and for the case of B n and F4, 
an extra prescription is required.) We then follow the (standard) procedure |25|]of defining 

Z { f\X) = tg(A) 4EE E Q % b) ^ A " W )' ( 8 ) 



N 

6=1 k=l v=l 



so that the (half-) integers 1-^ satisfy 



Zf\x^) = 1 ^. (9) 
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If we assume that Zj a \\) is monotonic, then the range of the integers is set by taking the 
difference of the limiting values, 2/j°'(±oo). For a given choice of the phase factor Q^(= 
e (2nsi/L) wnere ground state lies), the first striking consequence of these assumptions 
is that for strings of length I, the range of integers Aij coincides with the total number 
of strings of that length, , corresponding to "no holes" as seen in the bulk calculations 
of |Q] and ||. This reinforces our understanding that all the statements made in the bulk 
using densities can be upgraded to finite lattices. This will be exploited further later. Since 
the /-strings do not contribute to the counting of states, we can eliminate Mj using the 
sum rule @) . The range of integers associated with the other strings are: 

A/j a) = NS^CJJ,, + Mf - £ ^(CgUiC^J jk M$\ j < /, (10) 

6=1 k=l 

where Cg is the Cartan matrix of the (simply-laced) Lie algebra Q and C^ 1 is the inverse 
Cartan matrix of This gives a combinatorial count of the number of states of the 

form described above, 

/A/ (a \ 

e n M <«> • <»> 

{M («0} 3,a KIV1 J 7 



We shall use the equality of binomial coefficients 

(12) 



A )= A 
B \A-B 



to define a new variable 



N^ a) = Aij a) -Mj a) , (13) 

which counts the number of holes. We shall have some more to say about this shortly. 

To make the idea of upgrading bulk quantities to finite lattices concrete, we define 
densities for strings and holes |]ll]]by requiring their integrals over all A to be exactly equal 
to the number of strings and holes divided by N, the size of the system. This is the precise 
statement that rules the interpolation between the bulk and finite lattices. In the notation 
of i, 

p( a) (0) = Mj a) /iV, ^ a) (0) = ivj a) /iV, (14) 

where /(0) denotes the Fourier transform of / for zero argument. In the thermodynamic 
limit, because solving the integral equations involves taking derivatives on the densities, 
one can choose the branches of the logarithms (of the Bethe equations) freely. What we 
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see is that there exists a choice of branch on the finite lattice that makes the translation 
of statements true in the bulk to finite lattices possible. 

Instead of setting up the equations for the thermodynamics here, it is more convenient 
to directly quote the result from ref. ||. We perform the trivial steps of setting x (the 



Fourier transform variable) to zero, and applying the definition (|l4])to eqn. (2.16) of ref. 
§to get 

ivf = N5 ap [C- A ] vi J 3S -EE^f, (15) 

6=1 k=l 

where 

K tb = (a a \a b ){Toam(t b j,t a k) - }. (16) 



jk _ , .. | , f...;.. ,j_ ; j_ j^ - 

ab 



/ 

We can also solve for M-" ; in terms of iV w to get 



=NJ2 {K-^[C- A ] pi J ka -ED^W' ( 17 ) 

fc=l 6=1 fc=l 

where K~ x is defined by 

r t c Z-l 

£ E eOS^S* = ^ifc- (is) 

c=l m=l 

In the above, the index t a Z never shows up because there are never any holes in this sector, 
and the sum rule (|2])is used to eliminate the dependence on M^J . Note that for the simply- 
laced cases, K factorizes into a "level" and a "rank" piece, and (p~5|) reduces to fllCf) using 
©■ 

We can thus upgrade (|ll|)to the general case of all untwisted affine Lie algebras, where 
the number of states is of the form 



M (a) N (a) M (a) ^(o) 

J J \ _ I J J 

Mj a) J \ ivj a) 



(19) 



Once the counting problem is set up, it is possible to extract information about the ground 
and excited states by relying on the (tacit) assumption that the physical quantities of in- 
terest are "diagonalized in the string basis," without even specifying the dependence of 
the energy of the system on the roots of the Bethe equations! If the one-dimensional 
hamiltonians are derived from the transfer matrix of a 2-dimensional classical statistical 
mechanical model with positive Boltzmann weights, then the Perron-Frobenius theorem 
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forces the ground state to be unique. For the RSOS models, we are fortunate in having a 
situation where both the largest and smallest eigenvalue of the 1-dimensional hamiltonian 
correspond to regimes in the 2-dimensional model where the Boltzmann weights are posi- 
tive. Therefore the ground states for the hamiltonians with either a positive or a negative 
overall sign are unique. 

One candidate is, of course, the sea of t a l strings, because it has no holes. For the 
other candidate, we have to solve for the equation = which implies that the number 
of strings in that state is 

t p i-i 

M^=Nj2(K-^ p [C^ pi J ks . (20) 
fc=i 

Since the ground state must, by definition, be stable under an arbitrarily large but 
finite number of excitations, we can see that the natural variables to discuss the excitations 
are the Ms for the ground state of t a l strings, and the iVs for the case where the ground 
state is charecterized by the solution to (pCf ) . The momentum eigenvalue is written as a sum 
of the left hand sides of (f|) @ || (as in Bethe's solution of the isotropic spin 1/2 spin chain 
[Ij]). Thus, the sum of the integers for any set of roots gives the total momentum (with the 
appropriate ^ taken into account) of the state that these roots correspond to. Therefore 
we see that for every term in AI^ = + for which the term proportional to N 
does not have support, the contribution of these strings (or holes) to the momentum of the 
eigenstate at order one (N ) must necessarily vanish. For a system with no mass gap the 
energy contribution must consequently be zero to order one. However, they could (and do) 
contribute to the spectrum at order 1/N. These have been termed "ghost" excitations in 
||13|| . The strings (and holes) whose ranges are macroscopic (i.e. have a term proportional 
to N) constitute the order one excitation spectrum, and the coefficient of N encodes the 
Brillouin zone scheme of these excitations. 

To make these discussions more concrete, let us specialize to the simply-laced case 
(all t a = 1) where K^ b = (Cg) a b[C^_ i ]jk- Recall that for one of the two regimes (e = +1) 
in Qand 0, the state filled with / strings is the ground state, and all the other strings 
which have a macroscopic range (i.e., for which the Kronecker delta in the term on the 
r.h.s. of ([TOD proportional to N has support) have non-zero energy. As an example, the 
3-state Potts model is in the same universality class as (and in fact can be obtained by 
orbifolding) the RSOS model corresponding to SU(2), with I = 4. We can read off the 
Brillouin zone scheme of ||from eq. (|T0|). (Technical aside: also of interest is the direct 
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correspondence between the "zeroes of Q" as studied here, and the "zeroes of T" as studied 
in ||, ||. For I = 4, 2Mi corresponds to the number of 2-strings, 2M-i to the number 
of ns, each dressed by 2 plus roots, 2M3 to the number of —2-strings each dressed by 4 
plus roots and 2M4 to the (even) number of minus roots each dressed by 3 pluses, in the 
terminology of ||. The counting works out analogously, as does the correspondence with 
the ground state structure and excitation spectrum. It would thus be of interest to study 
this in more detail.) 

For the other choice of ground state, which corresponds to the case e = — 1 in 0, 
we find, from eqn. Q2"U|), that the state with only Ms^ non-zero has no holes, and is the 
vacuum of the theory, and these are the only type of strings that contribute to the energy 
and momentum at order 1 (in N). In the relevant hole variables, the integer ranges look 
like 

Alf =N8 sj [C^]+N^ -tWU^jrf. (21) 



and makes the last remark transparent. 

At this point we make a note that following the correspondence of [|34[], if the string 
hypothesis gives an exhaustive count of all the states in this (ground state) sector, these 
sums of products of binomial coefficients are precisely the multiplicities of the conformal 
blocks made up of tensoring iV copies of su> p representations and projecting on the singlet, 
and can be computed using Verlinde's formula [ Z7|j . The inhomogeneous chain will give 
the analogous count for tensoring arbitrary representations. This has also been observed 
in |32[]for su(2). In the simplest case of the fundamental representations of su(2), these 



numbers are indeed the number of paths leading to the singlet on a truncated Bratelli 
diagram as in eq. (1.21) of [[37]] (see also [p8[ ). The binomial sums can be performed 
in terms of generating functions, and the answer for su(2) (following the Appendix of 
|2~4| ) is expressed in terms of contour integrals over Chebyshev polynomials, suggesting 



a possible connection with the fusion rings of [4~3|. Not surprisingly, the pole terms that 
contribute are precisely those algebraic numbers that satisfy the infinite spectral parameter 
thermodynamic Bethe ansatz equations that determine the arguments of the dilogarithms 
that give the central charge |H] || [jDJ . For the higher rank algebras, we expect to encounter 
the group characters (following P3|| ). However, this analysis is still incomplete and we hope 
to report on this in the near future, where we shall also present some details on Yangian 
characters. (It is easy to see that the formulas in the appendix of ||can be obtained 
from our combinatoric formulas by setting I to infinity, which is the "classical" limit of 
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the "quantum group." With these explicit formulas the truncation condition in ||can 
presumably be verified.) 

We are now in a position to construct the branching functions for affine Lie algebras. It 
is seen in numerical studies (which provided the background for ref. |J), that within each 
class of states with a paricular root content, the sum of the absolute values of the integers 
give a good estimate of the (approximate) degeneracy of the levels, i.e., those states with 
the same value for the sum of the absolute value of the integers had energy eigenvalues that 
were almost equal. This is reminiscent of the conformal field theory definitions of energy 
and momentum being the sum and differences of the eigenvalues of Lq and Lq. One could 
then presume that this degeneracy would become exact in the thermodynamic limit, and 
counting degenerate states could be achieved by keeping track of these integers. There 
are important variations on this observation, however, in that the set of strings that do 
not have macroscopic ranges actually contribute both a positive and a negative amount in 
the energy eigenvalue at order 1/N. It is because these integers (actually, the momenta 
that correspond to them) characterize the "fermionic" (distinct integers!) quasiparticle 
spectrum, the resulting modular forms constructed by (/-counting or, alternatively, building 
up the partition function in a particular sector, are called "fermionic" in f[2l , [|13| , |16| . 



We shall distribute these integers with powers of q by using the the symmetric (under 
q — > g _1 ) (/-numbers (as in p6fl ) 

(g n/2 _ -n/2) 
W = ( g l/2_ g -l/2) ' U Z ' 

n 

W! = IIW' (22) 

i=l 



and construct the objects 



B) a {B}\ 



{M {a) } 3> a 3 



A/j a) 

{M (a) } 3> a 3 



(23) 



*= e n 



In order to focus on either chiral half of the spectrum, it is more instructive to write 
these objects in terms of g-binomial coefficients, which are polynomials in q and are defined 

as 



where 



(q;q). 



A 



(q;q)A-B(q;q)B 

A 
3=1 



(24) 



(25) 



(non-zero only for integers, A and B, with < A < B). These are then appropriately 
symmetrized under q — > q~ x by dividing by one-half of the degree of the polynomial. We 
thus have 



\/ Y, Il l 

{M (a) } 3, a 



±M (a) N (a) 



Xh 



Err -hM {a) N {a) 
[[q 2 ' 3 



{M (a) } 3>a 



Mj a) + ivj a) 
Mj a) 

M (a) + N (a) 



(26) 



N 3 



(a) 



where we write \ v only in the variables and Xh only in the ivj a ^ variables, using 

(Ifopand (|17|) . Since large values of A correspond to small momentum contributions and 
therfore, low energy states, and the integers have a monotonic dependence on A, the 
largest integers correspond to the lowest energy, and we assign zero momentum to the 
(half-) integer N/2 that lies in the power of q multiplying the g-binomial coefficients. 
Now the Xp an d Xh look like 



Xh 



hM-K-M 



n 



3,a 



i 1 

{M ( . a) } 

q-^ N < K ^> N x 



US K js 4- M (a) 4- V r V ibZ_1 K jk M { 

2 ap J\ ap -|- 1V1- -rl^b=ll^k=l IV ab 1V1 k 



M) 



(a) 



{N (a) } 



(27) 



x 



n 



3, a 



t b i-i, 



-l\jk AT (b) 



r(a) n 



N 



(a) 



with the condition that 



lim 



N 
m 



I 



(?; q)i 



(28) 



These are the branching functions (in the vacuum sector) corresponding to the con- 
formal field theories constructed as cosets ( a la GKO) p9] . The list of branching functions 



these (/-series expansions correspond to are given in [EJ. The central charges which can 



be obtained by taking the q — > 1 limit as in ||20|| , [lTjand [|T5|]are also listed there. One 
can read off a large number of non-trivial equivalences between models based on different 
groups and at various levels of fusion and roots of unity (level-rank dualities |41|] j42| ) . To 
come back the 3-state Potts again, for example, we see that the G% level 1 model for 
p = 2 (arbitrary s) gives the Lepowsky-Primc ||40|1 , ||14| , |l8|]forni for the c = 4/5 character. 
One could surely find lots of other interesting examples. 

It is amusing that the q — ^ 1 limit for the symmetric (/-binomial (the {}form) gives 
the total count of the states, while by pushing the focus onto the low-lying states (the [ ] 
form), one can extract the entropy of the low lying states in the same limit. 

In [PD| , there is a simple prescription for restricting to the type II [pJTf |39[[integrable 
representations which are relevant for these models - namely, the roots of the Bethe equa- 
tions are non-singular. They also propose that in order to go to a different sector one 
has to send a string to infinity. The only change in the above formulas (in the M lan- 
guage) would be to shift the Mj variable by Sji, where i refers to the particular length 
of the string that is sent off to infinity. This immediately gives the Lepowsky-Primc form 
| 4T[||111 |[18|for all the A± branching functions (after setting Q to su(2)). The restriction on 
the summation variables reflect the total number of strings in each sector. We conjecture 
that the same procedure for the higher rank cases would give all the branching functions 
of the generalized parafermionic theories. The transformation of this operation to the hole 
language would give the branching functions in the other regime, but the exact principle 
behind the restrictions on the integers seems unclear so far. 



It is intriguing that in ref ||19|1 , the authors produced their character formulas by 
considering the branches of the dilogarithm, whereas we seem to do so with the ordinary 
logarithm(ic form for the Bethe equations). Further relations between these two approaches 
might lead to a greater understanding regarding the solutions of the Bethe equations. 

For the elliptic case considered in (§], the essential steps of setting up the g-series 
representation can easily be carried through. However, the presence of a mass gap in 
the system distinguishes between the contributions of the macroscopic and microscopic 
excitations to the partition function. While, for example, this observation does indeed put 
the Eg based character formulas and the structure of the integrable scaling field theory of 
the Ising model in a magnetic field in the same framework, a better understanding of this 
(/-series away from criticality is lacking. We hope to understand how the mass spectrum 
may be extracted by upgrading the combinatoric arguments to the massive case, as well 
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as some understanding of how this character construction is related to Baxter's corner 
transfer matrix method. It should also be noted that the thermodynamic Bethe ansatz 
equations may be set up for models which do not necessarily possess lattice Bethe ansatz 
integrability, and (/-series modular forms may thus be associated with any such model 
possessing a factorizable S'-matrix. 

These combinatoric formulas have deep relations with several objects of mathematical 



interest. In |35], bijective correspondences between the combinatorics of Young tableau 
and the counting of Bethe ansatz states in su(n) invariant systems have been proposed. In 
particular, they could identify an appropriate Lascoux-Schiitzenberger charge of tableaux 
(see [[36|] [[35] for definitions) in terms of the integer prescriptions to string solutions (what 
they called rigged configurations) to get expressions for g-Kostka polynomials (||36|1). The 
results presented here indicate that for restricted tableaux (with no more than / columns) , 
this counting procedure relates the g-Kostka polynomials of these infinite (for the filled 
fermi sea) tableaux with modular forms associated with affine Lie algebras. The connection 
with the combinatorics of crystal bases (an explicit construction has again been carried 
out via a Fock space realization on extended Young tableaux j45|) should provide insight 
into the relationship between the counting of solutions to the Bethe equations and the 
corner transfer matrix calculations mentioned above. There are, of course, several open 
questions we have mentioned, and several directions of continuing research. It would be 
interesting to see if the counting scheme is actually of some universal form and whether 
the classification of states of other models may be brought into this form. Two interesting 
cases are the Z n parafermions studied in [|23[] (through the "zeroes of T" of the Fateev- 
Zamolodchikov model |2^]) and the model of ||21|| . The character formulas (at least for the 
Z n models can, of course, be obtained by level-rank duality |15| . A study of the model 



in ref. []21[ should yield interesting results. The question being posed is whether there 
is information beyond what is implicit in the formulas for the branching functions that 
distinguishes between different models. The outstanding question is, of course, the validity 
of the string hypohesis - namely, do the solutions of ([!]) actually fall into these classes? If 
not, how does the counting procedure rearrange the non-stringy effects to produce such 
clean results? 
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